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Supersymmetry (SUSY) is one of the most interesting and comprehensively studied models for new physics 
beyond the Standard Model. If SUSY exists in nature the Large Hadron Collider will provide excellent oppor- 
tunities to search for SUSY. SUSY discovery strategies of the ATLAS and CMS experiments are presented with 
a focus on early data. SUSY mass measurement techniques and determination of SUSY model parameters are 
also demonstrated. 



1. Introduction 

Supersymmetry is a theoretically attractive model 
for extention of the Standard Model (SM) It pro- 
poses that all the SM particles have SUSY partners 
with a spin difference of ±1/2. Among its motivations; 
SUSY solves the fine-tuning problem (Higgs mass sta- 
bilization against loop corrections), SUSY modifies 
running of the SM gauge couplings just enough to 
give the grand unification at a single scale and SUSY 
also offers a candidate for dark matter through R- 
parity conservation. Since SUSY partners of the SM 
particles have not been observed at the same mass 
scale SUSY must be a broken symmetry at low en- 
ergy. The Minimal Supersymmetric Standard Model 
(MSSM) brings more than a hundred free parameters 
into the theory, thus searching for SUSY is a very 
challenging task. 

One of the main objectives for building the Large 
Hadron Collider (LHC) is to discover SUSY. The AT- 
LAS (A Toroidal LHC Apparatus) and CMS (Com- 
pact Muon Solenoid) are two general purpose detec- 
tors at the LHC designed to search for SUSY particles 
(and for the Higgs boson). Both experiments have 
successfully seen the first beam events in September 
2008 and have been commissioning with the cosmic 
rays since then. The LHC is scheduled to start up in 
November 2009. 

The ATLAS and CMS experiments put enormous 
efforts recently to get ready for the analysis of the 
first LHC data. ATLAS has written a collection of 
notes published as a CERN book Q based on the 
results from the Computer System Commissioning 
(CSC) exercise. CMS has written a Technical De- 
sign Report II P. Both have dedicated chapters to 
SUSY covering a broad spectrum in SUSY analyses 
with 14 TeV center of mass energy. Post-CSC stud- 
1 ISJiJS 9 and post-TDR studies 
111 . [l2l 1 131 ] have followed for both 14 
TeV and 10 TeV center of mass energies resulting in 
a variety of public notes. This review article is based 
on these recent studies. After a brief mention of the 
SUSY search strategies at the LHC, examples of the 
data-driven methods will be given for estimating the 
SM backgrounds to SUSY searches. Then searches 
with the inclusive missing transverse energy (E™ lss ) 
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signatures and an exclusive search without explicity 
use of E!f lss will be summarized. Examples from spe- 
cial SUSY signatures of the non-standard SUSY mod- 
els will be presented together with the discovery reach 
of some SUSY models. The article will conclude af- 
ter a discussion on the SUSY mass measurements and 
determination of SUSY model parameters upon a pos- 
sible SUSY discovery. 



2. SUSY Search Strategies at the LHC 



The SUSY breaking mechanism determines the phe- 
nomenology and search strategy in the collider ex- 
periments. Many possible mechanisms have been 
proposed; mSUGRA 111, GMSB [TBLAMSB [lj, 
SO(10) SUSY GUTs [lj Split-SUSYpj] and others. 
The search strategies can be discussed in two classes 
of models: 

1. Standard SUSY models with R-parity conserva- 
tion. In this class SUSY is searched for as an 
evidence of excess above the SM background in 
the following channels; 0, 1, 2 leptons+> 2, 3, 
4 jets+£™ ss , 7 -r-jets+ J Bf liss , r or b-\ets+E^ iss , 
multi-leptons, multi-jets without E^ lss and oth- 
ers. 

2. Non-standard SUSY models with special signa- 
tures; displaced vertices, stopped gluinos in the 
calorimeter, non-pointing photons, long-lived 
stable massive particles and others. 

Understanding the SM backgrounds play an es- 
sential role in SUSY searches. Main sources of the 
SM backgrounds come from the following processes; 
tt+jets, TU+jets, Z+jets, QCD jets and diboson (ZZ, 
WW, WZ) production. The data-driven estimation 
of the SM backgrounds, will be an emphasis of early 
analyses. Understanding the fake E™ lss and tails in 
the E™ lss distributions as well as the reconstruction 
of the non-standard signatures will be a critical task 
in the search for SUSY. 
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3. Estimation of the Standard Model 
Backgrounds to SUSY Searches 

Several Monte Carlo-driven and data-driven estima- 
tion techniques have been developed to estimate the 
SM backgrounds. Some of the data-driven estimation 
techniques are listed as: 

• Z — > vv background in 0-lepton mode [l3j . 

• QCD background in 0-lepton mode SUSY search 
using the jet smearing method [2|. 

• QCD background in 1-lepton mode SUSY search 
using the lepton isolation method Q ■ 

• tt background in 0-lepton mode using the re- 
placement technique [g]. 

• tt background in 1-lepton mode SUSY search 
using the Topbox method [2j. 

• Combined background in 1-lepton mode SUSY 
search using mx method 0] (further develop- 
ment with combined fit) and Tiles method [f|. 

• Others. 

As an example two of the data-driven approaches 
are summarized below. 

3.1. Data-driven estimation of the QCD 
background in the hadronic channel 

The biggest background estimation challenge for 
SUSY searches in the jets+£ , ™ lss channels will be esti- 
mation of the QCD jet background. This background 
involves two sources of E^ lss : 

1. Fake E™ lss ; non-Gaussian tails to the detector 
jet response function arising from the dead ma- 
terial, jet punch-through, pile-up and other ef- 
fects. It can be suppressed by applying the 
jet-_E™ ISS azimuthal angle cuts, calorimeter and 
tracking cuts and cosmic background rejection 
cuts. 

2. Real £J™ ISS ; resulting from the non-interacting 
particles such as neutrinos or the Lighest Su- 
persymmetric Particle (LSP). 

Two approaches are taken to estimate the remain- 
ing QCD jet backgrounds; Monte Carlo-driven esti- 
mates and data-driven estimates. The Monte Carlo 
based estimates are subject to large systematics ef- 
fects arising from the proton parton distributions, un- 
derlying event uncertainties, jet energy scale uncer- 
tainty and an uncertainty in modeling the QCD jet 
physics with Monte Carlo generators (PYTHIA ver- 
sus ALPGEN). There will also be an uncertainty from 



the luminosity measurement. In addition, the detec- 
tor simulation uncertainties (imperfect description of 
the response of the detector to jets) and the statisti- 
cal uncertainties (large QCD cross section rendering 
production of full simulation samples unfeasible) will 
contribute. As for the data-driven estimates they will 
be crucial in early phase of data taking as they are 
less prone to input systematics. 

A method from an ATLAS study @ is described 
here for estimating the QCD jet background in the 
0-lepton mode SUSY search. The method makes 
use of smearing the jet transverse momenta in low 
E™ lss QCD jet data with a data measured jet re- 
sponse function R (ratio of the measured jet pt to 
the true jet pt)- The Gaussian part of the jet re- 
sponse function is measured with the balance of the 
7+jet events. The non-Gaussian part of the jet re- 
sponse function is measured from the 'Mercedes' type 
of configuration (E T nss is parallel or anti-parallel to 
the pt of one of the three jets selected in the event). 
The jet response function is then used to smear the 
jet pt in multijet events with low £?™ ss in order 
to estimate the i?™ lss distribution of QCD multijet 
events. Figure [1] shows the _E™ 1SS distribution of the 




Er [GeV] 

Figure 1: The £"™ lss distribution of signal and background 
for an integrated luminosity of 23.8 pb" 1 at 14 TeV. 
The QCD background is estimated from the jet smear- 
ing method in 0-lepton SUSY search as explained in the 
text. The signal point is the ATLAS mSUGRA benchmark 
point SU3. 

full simulation data and the data-driven estimated 
for an integrated luminosity of 23.8 pb _1 at 14 TeV. 
The E^ iss distribution of the SUSY signal (ATLAS 
mSUGRA benchmark point, SU3, with parameters 
m = 100 GeV, m 1/2 = 300 GeV, A = -300, 
tan(/3) = 6, sgn(/i) = +) and the summed background 
from Z — > yp+]ets, W — ► Z^+jets, tt+jets are also 
shown. A good agreement is seen between the es- 
timated and the full simulation data. Event counts 
yield 2.36 ± 0.09(stat) ± 1.44(syst) for the estimated 
versus 1 for the observed. The uncertainty on the 
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estimate is then 60% for 23.8 pb^" 1 . The same un- 
certainty is assumed for 1 fb _1 of data. Given the 
difficulty of obtaining accurate estimates of the QCD 
background, results from a number of independent 
techniques should be compared to obtain a robust 
QCD background estimate. 

3.2. Data-driven background estimates 
for SUSY diphoton search 

The GMSB model with two high pt photons and 
large E™ lss in the final state is among the promis- 
ing models for early SUSY searches. A method 
from a CMS study [llj is described here to predict 
the E^ lss distribution in a diphoton sample from 
the SM processes. Seeing an excess of events at 
high E^ lss would indicate a signal for new physics. 
The physics background is small and comes from 
the Z"/"f — ► i/wyy and VU77 — * Iwyy SM processes. 
The instrumental background comes from three ma- 
jor sources: 

1. QCD background; results from 7-jet misidenti- 
fication in QCD events with no real E™ lss such 
as photon plus jets and multijet production. 

2. Electroweak background; results from 7-e 
misidentification in events with real E™ lss from 
Wj and W+jet production. 

3. Non-beam background; results from high energy 
muons from cosmic rays or beam-halo. 

It has been demonstrated that the QCD and elec- 
troweak backgrounds are under control, and the tech- 
niques explored will effectively eliminate the non- 
beam background. A CMS benchmark point, GM1C, 
is used in this study with parameters A=100 TeV, 
M m =200 TeV, N 5 =l, C G =1, tan(/3) = 15, sgn(» = 
+. Figure [H shows the i?™ lss distribution from the 
background closure test (using Z — > e + e~ events to 
describe the QCD background) in absence of SUSY 
signal for an integrated luminosity of 100 pb^ 1 at 10 
TeV. As for comparison with the Monte Carlo truth, 
event counts yield (2.78 ± 0.24) for the Monte Carlo 
77 data and (2.69 ± 0.66) for the estimated back- 
ground in absence of SUSY signal, and (17.5 ± 0.26) 
and (2.99 ± 0.68) respectively in presence of SUSY 
signal. A good agreement is seen between the data- 
driven estimated and the predicted background. The 
data-driven strategy demonstrated in this study can 
be used at the start-up of the LHC. 

4. SUSY Searches with Inclusive 
E™ ss Signatures 

SUSY phenomenology varies significantly through 
the SUSY parameter space, therefore SUSY is 
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Figure 2: The E™ iBB distribution of the Monte Carlo data 
and the data-driven estimated from a background closure 
test in absence of SUSY signal. An integrated luminosity 
of 100 pb" 1 at 10 TeV is considered. 

searched inclusively based on very general signatures. 
It is expected that gluinos and squarks (heaviest 
SUSY particles) will be produced in the initial proton- 
proton interaction at the LHC and then their cascade 
decays will result in events with high px jets and lep- 
tons in the final state. In addition, if the R-parity 
is conserved the LSP would be stable and undetected 
leading to events with large E™ lss . Thus typical in- 
clusive SUSY searches cover a broad range of search 
modes depending on the presence of leptons (elec- 
tron/muon), r lepton or b quarks in the final state; 
the 0-lepton mode, 1-lepton mode, 2-lepton mode, 3- 
lepton mode, tau-mode and 6-jet mode and others. 
Examples of these search modes are given in the fol- 
lowing. A review of the SUSY searches with inclusive 
E^ lss signatures from an ATLAS study can be found 
in more detail in this volume [l9[ . 

4.1. The 0-lepton and 1-lepton modes 

The 0-lepton mode (jets+i?™ ISS ) provides the least 
model dependent search. Events with at least four 
jets or three jets or two jets can be chosen. The 
main background in this mode is the QCD multi- 
jet background (Elf lss is produced either by fluctu- 
ations in jet energy measurement or by B hadron de- 
cays to real neutrino). A strong event selection cuts 
(requirement of high multiplicity jets, high cuts 
on jets and a cut on the azimuthal angle between 
each of the first three jets and the E™ lss ) reduces this 
background significantly. The remaining background 
comes from rt+jets, VF+jets, Z+jets and diboson pro- 
duction. The global event variable, M c g, is used in 
discriminating SUSY from the SM events. M e g is de- 
fined as the scalar sum of the pt of the leading jets 
and the E™ lss and the pt of the leptons (if present) 
considered in the analysis. M e s is also used to quan- 
tify the SUSY mass scale. From an ATLAS study 
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the M eS distribution of the signal (ATLAS mSUGRA 
benchmark point, SU4, with parameters mo = 200 
GeV, mi/2 = 160 GeV, A = -400, tan(/3) = 10, 
sgn(/i) = +) and the total background are shown in 
Figure [3] for an integrated luminosity of 200 pb _1 at 
10 TeV. An excess of events is clearly observed in this 
mode. 
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Figure 3: E™' 1BB distribution in the 0-lepton mode SUSY 
search for an integrated luminosity of 200 pb _1 at 10 
TeV. The open circles shows the signal (ATLAS mSUGRA 
benchmark point SU4) plus the total SM background while 
the different SM backgrounds are indicated in the legend. 

The 1-lepton mode (jets+one lepton+£?™ lss ) is of- 
ten considered to be the golden mode for early SUSY 
searches. Requirement of the additional lepton (elec- 
tron or muon) strongly reduces the QCD multijet 
background. The event selection cuts are similiar to 
that of the zero-lepton mode but includes a cut on 
the transverse mass of the lepton and the _E™ ISS to 
suppress the tt+jets and PF+jets backgrounds. Fig- 
ure 2] shows the E™ lss distribution for the same AT- 
LAS study as in the 1-lepton mode. Although the 
statistics is reduced, the background is dominated by 
tt+jets and W+jets which are expected to be bet- 
ter understood than the QCD background. Thus 1- 
lepton mode is more robust against uncertainities in 
the background estimations. 



4.2. The 2-lepton and 3-lepton modes 



The 2-lepton mode with jets+two leptons (opposite- 
sign or same-sign) +i?™ lss in the final state can have 
high discovery potential in early SUSY searches. The 
opposite-sign (OS) dileptons arising from the SUSY 
decays (x° — ► Xi^ ) have the same-flavor to avoid 
fj, — > e*y decay. The dileptons coming from different 
decay chains (SUSY combinatorial and the SM back- 
grounds) can be of the same-flavor (SF) or opposite- 
flavor (OF) with the same probability. Thus observing 
an excess of OSSF dilepton events over OSOF events 
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Figure 4: _B™ 1SS distribution in the 1-lepton mode SUSY 
search for the same ATLAS study as in Figure El 



would indicate a clear evidence for new physics. Fig- 
ure [5] shows the invariant mass of OSSF dileptons 
from a CMS study [To| for an integrated luminosity 
of 200 pb" 1 at 10 TeV. The CMS mSUGRA bench- 



mark point, LMO (to = 200 GeV, 



'1/2 



160 GeV, 



^4o = —400, tan(/3) = 10, sgn(/i) = +) is considered. 
The red histogram shows the contribution from SUSY 
events. The dominant background is found to be 
ti+jets. A significant excess of OSSF dileptons events 
over OSOF events (shown by the black solid line) is 
seen with the first 200 pb" 1 of the LHC data. The 
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Figure 5: The invariant mass of OSSF dileptons for an 
integrated luminosity of 200 pb _1 at 10 TeV. The red 
histogram shows the SUSY signal while the different SM 
backgrounds are indicated in the legend. The black solid 
line respresent the background from the OSOF pairs. The 
black dots shows the distribution when exactly 200 pb _1 
of Monte Carlo events are used, namely without scaling 
down from a larger sample. 

same-sign dileptons could provide an almost back- 
ground free channel for SUSY searches. Events with 
two prompt, isolated and same-sign dileptons are rare 
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in the SM while in SUSY the gluino is a Majorana par- 
ticle so events with same-sign dilcptons can be com- 
mon. Though the event rates will be very small (from 
10 to 100 events/ftr 1 ) a recent ATLAS study || dis- 
cusses that gluinos and squarks with masses up to 
500 GeV can be discovered with a 3a signal signifi- 
cance for an integrated luminosity of 200 pb _1 at 10 
TeV. 

The 3-lepton mode (three leptons+£'™ lss ) to search 
for signals from direct gaugino (such as \\ an d Xi) 
pair production has perhaps been a best search mode 
at the Tevatron [2(|. The trilepton signal can come 
from leptonic decay of pairs of heavy gauginos (pro- 
duced directly or in the decay of heavier partner par- 
ticles) through real or virtual W^, Z° or I to leptons 
and a pair of LSPs. A general approach has been to 
consider two analyses; first with the selection of the 3- 
leptons+jet events (no cut on E™ ss ) and second with 
the selection of the 3-leptons+i?™ ISS events. In an AT- 
LAS study [H, the 3-leptons+jet analysis is found to 
be more sensitive than the 3-leptons + E^ lss one for 
the mSUGRA benchmark points studied for an inte- 
grated luminosity of 1 fb _1 at 14 TeV. 




Figure 6: The M e e distribution of the b-jet mode (two 
b-jets analysis) SUSY search for an integrated luminos- 
ity of 1 fb" 1 at 14 TeV. The signal point is the ATLAS 
mSUGRA benchmark point SU6. The contributions from 
the signal and the SM backgrounds are indicated in the 
legend. 



5. An Exclusive Search without E^" 



4.3. The tau- and 6-jet modes 

As for the tau- and &-jet modes, they can be dom- 
inant in SUSY models with large tan(/3), thus they 
can help determining the underlying model param- 
eters. The requirement of a reconstructed r in the 
final state eliminates the QCD background and the 
remaining background is dominated by tt+jets and 
JU+jets. Though the lower efficiency and purity in r 
reconstruction, r decays provide complimentary info 
on the nature of new physics; f LSP, e/ li/t universal- 
ity. 

In mSUGRA models with large tan(/3) the b and 
t are lighter than the first and second generation 
squarks. In this case b production is enhanced lead- 
ing to decays with &-jets in the final state. A similiar 
analysis is performed as in the 0-lepton mode with 
the additional selection of at least two or three 6-jets. 
The 6-tagging performance at high px is an important 
issue for identifying &-jets. A recent ATLAS study [H 
uses a ^-tagging efficiency of ^60% and a rejection 
of ~100 (~10) against light quark (c-quark) jets and 
shows that requiring 6-jets reduces the QCD multi- 
jets, W+jets and Z+jets backgrounds leaving the al- 
most only background from ti+jets. Figure [5] shows 
the M c e distribution of the mSUGRA signal point 
SU6 (m = 320 GeV, m 1/2 = 375 GeV, A = 0, 
tan(/3) = 50, sgn(^i) = +) together with the SM back- 
grounds for an integrated luminosity of 1 fb -1 at 14 
TeV. A clear excess is seen for the signal events above 
the SM backgrounds. 



A new approach to SUSY searches with dijet events 
has recently been proposed [2l| . It focuses on a SUSY 
parameter space where squarks are pair produced and 
both decay directly to a quark and the lightest neu- 
tralino, leading to dijet events with E™ ss in the final 
state. It exploits powerful discriminating variables to 
separate signal and background without making ex- 
plicit use of the E™ lss measurement, thus providing a 
very promising channel for early SUSY searches. A 
CMS study based on this approach is summarized 
in this section. The study has been carried out in 
the context of the mSUGRA model for an integrated 
luminosity of 100 pb" 1 at 10 TeV. 

The main backgrounds to the event topology are 
the QCD dijet events (i? T 1ISS is introduced through 
jet energy mismeasurements) and Z+jets events with 
Z decaying into two neutrinos. A kinematic vari- 
able, ax — E^/Mt, is defined as the ratio of the 
transverse energy of the second jet and the transverse 
mass of the two jets. It exploits the requirement of 
back-to-back jets of equal magnitude for QCD events, 
thus ax is exactly 0.5 for well measured QCD dijet 
events. Figure [7] shows the Qt distribution for di- 
jet events for signal and background. The mSUGRA 
benchmark points are labelled as LMO (mo = 200 
GeV, m 1/2 = 160 GeV, A = -400, tan(/3) = 10, 
sgn(^) = +) and LM1 (m = 60 GeV, m 1/2 = 250 
GeV, A = 0, tan(/3) = 10, sgn(^) = +). Events arc 
selected with ax > 0.55 cut which removes almost all 
QCD background. The remaining small background is 
from tt+jets, IU+jets and Z — > w+jets events. This 
study also extends the dijet analysis to higher jet mul- 
tiplicities where n jets (n — 2. ..6) are considered. It 
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is found that the discriminating power of «t against 
QCD background provides signal over background ra- 
tios of 4 to 8 for favorable SUSY benchmark points. 
The results are shown to be robust against system- 
atic uncertainties from the jet energy scale and jet 
direction. A data-driven approach has been taken to 
estimate both the QCD background and the remain- 
ing backgrounds. Even though the above Monte Carlo 
study suggests that the QCD background is negligible 
(after the ar cut), it has been verified by the data- 
driven estimates that it is indeed the case. The details 
of this study can be found in this volume [24] . 




Figure 7: The ot distribution for dijet events for an inte- 
grated luminosity of 100 pb" 1 at 10 TeV. The labels LMO 
and LM1 represent the CMS mSUGRA benchmark points 
as their parameters given in the text. The SM backgrounds 
are also indicated in the legend. 



6. Special SUSY Signatures 

A number of SUSY models predict specific signa- 
tures that may not be observed by the above tradi- 
tional SUSY searches, such as long-lived massive par- 
ticles or high photons. For these specific signatures 
the SM backgrounds are often small thus their detec- 
tion could provide unique constraints on the SUSY 
breaking mechanisms. A list of SUSY models giving 
rise to stable massive particles can be found in [23]. 
Examples are MSSM, GMSB, AMSB and Split-SUSY 
models. 

These particles live long enough to pass through 
the detector or decay in it. Such a particle can be 
a slepton (ti, 6r, flu in GMSB model), a gluino (in 
Split-SUSY model), a stop (gravitino LSP scenario 
of SUGRA models [24() or a lightest neutralino (in 
GMSB model). Their signatures will be different at 
the detector. If it is a stable charged slepton, its in- 
teraction with the detector will only be ionizations, 
its observed track will look like a muon, except for its 



higher energy deposition and longer time of flight than 
a muon. If it is a gluino or a stop, it is meta-stable, 
it will form bound states, so-called R-hadrons. Their 
signatures will be similiar to that of the stable slep- 
ton together with the appearance of high j?t tracks in 
the muon system with no matching track in the inner 
detector or the electric charge flipping between the in- 
ner detector and the muon system. If it is a long lived 
neutralino, it can travel a significant distance before 
decaying into a photon and a gravitino. Such pho- 
tons will not point back to the beam interaction point 
called as non-pointing photons. If the neutralino's life- 
time is not too long, events with two high photons 
(prompt photons, originate close to the beam inter- 
action point) plus large E™ lss (from gravitinos) are 
expected. 

As an example one study on the detection of the 
long lived gluinos is summarized in what follows. The 
Split-SUSY model suggests a large mass splitting be- 
tween the new scalars and new fermions, then gluinos 
can only decay through a highly virtual squark if R- 
parity is conserved. Gluinos can thus be long-lived 
and may well be stable on the typical experimen- 
tal timescales. When produced these gluinos will 
hadronize into bound states (R-hadrons), gg, gqq, 
gqqq. A recent CMS study (T^j discusses a search 
strategy to detect the decays of such R-hadrons. The 
charged R-hadrons will lose energy via ionization as 
they traverse the detector, this energy loss may be 
sufficient to stop a significant fraction of the particles 
inside the detector volume. These stopped R-hadrons 
may decay at times when there are no collisions (beam 
gaps) or no beam (interfill period) in the LHC. Ob- 
serving such decays would immediately indicate a new 
physics signal. 

A custom Monte Carlo simulation has been devel- 
oped to estimate the signal efficiency and a novel 
triggering strategy has been implemented to improve 
the signal sensitivity. The only physics background 
is from cosmic rays. There is also an instrumental 
background. The uncertainties on the background de- 
termination are statistical and they are found to be 
small. However, the search sensitivity to a particu- 
lat model (Split-SUSY) involves some significant sys- 
tematic uncertainties; theoretical uncertainties in the 
NLO calculation of the gluino-gluino production cross 
section, uncertainty on the stopping efficiency from 
the Monte Carlo simulation, uncertainty on modeling 
the exponential decay of the instantaneous luminos- 
ity. Figure [8] shows the signal significance that can 
be achieved after 30 days of running with an instanta- 
neous luminosity of 10 31 cm -2 at 10 TeV as a function 
of the gluino mass from a counting experiment. A 5er 
discovery should be possible in a matter of days for 
the models with large cross sections (~1 nb). 
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Figure 8: The signal significance achievable after 30 days 
of running with an instantaneous luminosity of 10 31 cm -2 
at 10 TeV as a function of the giuino mass from a counting 
experiment. The giuino lifetimes (from one micro second 
to one week) are indicated in the legend. 
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Figure 9: The 5a reach contours for the 0-lepton, 1-lepton 
and 2-lepton search modes for mSUGRA as a function of 
mo and m\n. The horizontal and curved grey lines indi- 
cate giuino and squark mass contours, respectively. The 
scanned points are marked with a dagger. An integrated 
luminosity of 200 pb _1 is assumed at 10 TeV. 



7. SUSY Discovery Reach 

In the previous sections examples of SUSY search 
strategies have been given based on the analyses of 
SUSY benchmark points with R-parity conservation. 
These points can not be representative of all possi- 
ble SUSY models to be found at the LHC. They have 
been chosen with the aim of studying a variety of sig- 
natures and developing analysis techniques that can 
be applied to much of the SUSY parameters space ac- 
cessible with early LHC data. In order to map the 
5ct discovery reach of a SUSY model, a scan over its 
parameters is performed. 

A recent ATLAS study focuses on scanning the 
parameter space of mSUGRA, pMSSM (phenomeno- 
logical MSSM with 19 free soft SUSY breaking pa- 
rameters) [25| and UED (Universial Extra Dimen- 
sions) [261 ] models considering an integrated luminos- 
ity of 200 pb" 1 at 10 TeV. Figure M shows the 5a 
discovery reach as a function of mo and TO1/2 param- 
eters from a scan in mSUGRA with other parame- 
ters set to tan(/3) = 10, A = and sgn(/i) = +. 
The signal Monte Carlo samples were produced using 
ATLFAST2 simulation (GEANT4 simulation of the 
inner detector, the muon system and fast simulation 
of the calorimeter). The top pair, single top and dibo- 
son samples were produced using GEANT4 simulation 
while all other background samples with ATLFAST2. 
The reach plot has been made by finding the optimal 
M c ff cut to maximize the signal significance in the 0- 
lcpton, 1-lcpton and 2-lepton search modes. This plot 
includes uncertainties on the SM background deter- 



mination which are re-estimated from the data-driven 
methods studied in an earlier work 0] for an inte- 
grated luminosity of 1 fb _1 at 14 TeV. It is taken 
as 50% for the overall background uncertainty in this 
study. As seen in the figure gluinos and squarks can be 
discovered with masses up to 750 GeV. The 0-lepton 
mode is found to be the best mode, however q-lepton 
mode is more robust against the QCD background. 
The reaches for the 2-lepton modes are less than the 
reaches for the 0-lepton and 1-lepton modes. 

As for the 5cr discovery potential of the tau-modc, 
it is found to be slightly worse than the 0-lepton and 
1-lepton modes due to the lower efficiency and purity 
in t reconstruction which is studied in [2] for 1 fb _1 
at 14 TeV. 

The 5cr discovery potential of the 6-jet mode is 
found to be comparable to that of the 1-lepton mode 
from a recent ATLAS study U for 1 fb" 1 at 14 TeV. 
It is a competitive search mode specially at higher Too 
values. 

Scans in the NUHM (Non Universal Higgs 
Model) H3 and GMSB SUSY parameter space have 
also been performed in an ATLAS study for 1 ftr 1 
at 14 TeV. The discovery reach of the NUHM model 
is found to be vitually identical to that of mSUGRA 
in the 0-lepton and 1-lepton modes. The 5<7 discov- 
ery reach of the GMSB model as a function of tan(/3) 
and A parameters can be seen in Figure HH1 with other 
parameters set to M m =500 TeV, Af 5 =5, Cq=1 and 
sgn(^i) = +. The reach for the 3-lepton analysis is 
significantly better than that of the 2-lepton analysis. 
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Figure 10: The 5a reach contours of the 2-lepton and 3- 
lepton analyses for the GMSB SUSY parameter scan as 
a function of tan(/3 and A. The vertical solid and dashed 
grey lines indicate gluino and squark masses respectively. 
An integrated luminosity of 1 fb _1 is assumed at 14 TeV. 



8. SUSY Mass Measurements 

The searches strategies to establish SUSY discovery 
have been outlined in the previous sections. As soon 
as a SUSY discovery can be made by these searches, 
the emphasis will move on to confirming that it is re- 
ally SUSY. The masses and spins of the newly discov- 
ered particles should be measured and the underlying 
model parameters should be determined. Since it is 
not possible to cover all of the allowed SUSY models, 
mSUGRA model is taken as an example to develop 
the mass measurement techniques that can be appli- 
cable to large variety of SUSY and SUSY- like models. 
The study presented here is based on an ATLAS study 
performed for integrated luminosity of 1 fb _1 at 14 
TeV. 

The mass measurement strategy is to exploit kine- 
matics of the long decay chains. In the R-parity con- 
serving mSUGRA models SUSY particles are pair pro- 
duced and their cascade decays typically have high 
Pt jets, leptons and large E™ lss (due to two invisible 
LSPs in every SUSY event) in the final states. No 
mass peaks can be reconstructed, however kinematic 
endpoints in the invariant mass distributions of dilcp- 
ton, dijet, leptons+jets and ditau can be measured. 
These measurements are then used in deriving rela- 
tions between SUSY masses. The theoretical positions 
of the endpoints are well known from the analytical 
expressions [28j . If sufficiently long decay chains can 
be isolated and enough endpoints can be measured 
then the masses of the SUSY particles involved in the 
decay can be obtained in a model-independent way. 
As an example from the decay chain: 



the masses of gx, Ir, \2 anc ^ Xi can be measured. 
This is the first decay chain likely to be reconstructed 
as it has the advantage of having charged leptons, 
E^ lss and hadronic jets in the final state which en- 
sures a large signal to background ratio. Also, the 
technique known as flavor subtraction allows to deter- 
mine both the SUSY combinatorial and the SM back- 
ground from the data itself with high accuracy. This 
technique exploits the fact that the signal contains 
two opposite-sign same-flavor (OSSF) leptons, while 
the background is due to pair of leptons coming from 
different decay chains, which can be of the same-flavor 
(SF) or opposite-flavor (OF) with the same probabil- 
ity. Thus the background cancels in the substraction 
N(e + e~)//3 + /?A(^+/i~) - N(e ± ^) where N indi- 
cates the respective number of events and (3 is the 
ratio of the electron and muon reconstruction efficien- 
cies (/3=0.86 for this study). 



8.1. Kinematic endpoints 

Starting with the two leptons (electron or muon) 
in the final state, the invariant mass distribution of 
dilepton exhibits a triangular shape with an endpoint 
at: 



edge 



X 



1 - 



\ 



1 - 



(2) 



(1) 



In Figure [TT1 the dilepton invariant mass distribution 
is shown after flavor subtraction applied for the AT- 
LAS benchmark point SU3. The line histogram shows 
the small SM background contribution (dominated by 
tt+jets), while the points are the sum of SM and 
SUSY contributions. A clear excess is seen together 
with a clear edge structure. The distribution is fitted 
with a triangle smeared with a Gaussian which is su- 
perimposed in the plot. Also, the expected (truth) po- 
sition of the endpoint is indicated by a vertical dashed 
line. The fitted endpoint is (99.7 ± 1.4 ± 0.3) GeV 
where the first error is statistical and the second is 
the systematic error on the lepton energy scale and 
the f3 parameter (10%). Compared to its truth value 
of 100.2 GeV (calculated from Eq.[2]) the endpoint can 
be measured with a precision of a few percent already 
with 1 fb _1 of integrated luminosity at 14 TeV for the 
model chosen. 

In order to determine the masses of all the SUSY 
particles involved in the decay chain, further kine- 
matic endpoints need to be measured from the mass 
distributions involving a jet: m^ e , mf^*, m^f*. 

and m £q( S h ig h)- The label thr indicates lower end- 
point of the distribution (a non-zero threshold value) 
whereas low/high indicate minimum/maximum of the 
two masses mi+ q and mi- q . Only the two leading jets 
are considered in event selection as it is not possible 
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Figure 11: The dilepton invariant mass distribution for 
the ATLAS benchmark point SU3 with an integrated lu- 
minosity of 1 fb _1 at 14 TeV. Different distributions and 
lines are explained in the text. 



to identify the quark from the (Jl decay. Therefore 
it is assumed that it hadronizes into one of the two 
highest px jets in the event. The fitted endpoints to- 
gether with their truth values can be seen in Table |TJ 
The first error is statistical, the second and third are 
the lepton energy scale uncertainty and the jet energy 
scale uncertainty, respectively. While the measured 
values are compatible with their truth values, there 
are large uncertainties on the measured values. This is 
due to not having clear edge structure (combinatorics 
from choosing the wrong jet in a true SUSY event) 
and thus trying to fit the tails (beyond the endpoint) 
with straight lines. 



Table I Measured endpoint values for the ATLAS 
mSUGRA benchmark point SU3 with an integrated lu- 
minosity of 1 fb _1 at 14 TeV. The truth values are also 
listed. 



Endpoint 


SU3 Measured (GeV) 


SU3 Truth (GeV) 


edge 


517 ± 30 ± 10 ± 13 


501 


thr 

m Uq 


265 ± 17 ± 15 ±7 


249 


"%(low) 


333 ± 6 ± 6 ± 8 


325 


edge 
m ^(high) 


445 ± 11 ± 11 ± 11 


418 



In a similiar way, the endpoint of the invariant mass 
distribution of ditau, m°^. ge , can be measured from 
the decay x° — * — * Xi T±1 ~ T - The fitted endpoint 
is (102 ± 17 stat ± 5.5 syst ± 7 po1 ) GeV for the ATLAS 
benchmark point SU3. The first error is statistical, 
the second is the systematic uncertainty dominated 
by the fitting procedure and the third is the uncer- 
tainty from the polarization effects. Its truth value is 
98.3 GeV. 

Also, from the decay qn — > XiQ the mass of the 
right-handed squark can be measured from the end- 
point of the stransverse mass distribution by assum- 



ing that the mass of Xi is known from the dilepton 
and leptons+jets endpoint measurements. The fit- 
ted endpoint is 590 ± 9(stat)±g 3 (syst) for the ATLAS 
benchmark point SU3. The systematic uncertainty 
accounts for the choice of the fit limits and the jet 
energy scale variations. It agrees well with its truth 
value of 611 GeV. 



8.2. Extraction of SUSY masses from 
endpoint measurements 

From a combination of experimentally measured 
endpoints as above the SUSY particle masses can be 
extracted. By using the five endpoint measurements, 

edge edge thr edge edge n r 

m U ' m Uq I m U r "Vlow)' "Vhigh)' f ° Ur U1> 

known SUSY particle masses involved in the decay, 



5l, In, x\i Xi can be solved. A numerical \ 2 mini- 
mization based on the MINUIT package is performed. 
The masses resulting from the \ 2 fit can be seen in 
Table (Til for the SU3 point. The first error is sta- 
tistical, the second accounts for the jet energy scale 
uncertainty. An anticorrelation with the jet energy 
scale variation is indicated by ±~ sign. The recon- 
structed masses are found to be highly correlated with 
Xi which is not well determined by the endpoint tech- 
niques. Therefore the precision on the absolute mass 
values is rather moderate. However, the mass differ- 
ences are better measured than the absolute masses 
as expected. 



Table II Reconstructed SUSY particle masses and mass 
differences together with their truth values for the SU3 
benchmark point with 1 fb _1 at 14 TeV. 



Observable 


SU3 Reconstructed (GeV) 


SU3 Truth (GeV) 


m~ 


88 ± 60 T 2 


118 


m~ 


189 ± 60 =F 2 


219 


m~ 

1L 


614 ±91 ±11 


634 


rrtr- 
in 


122±61=f2 


155 




100.6 ± 1.9 ± 0.0 


100.7 


m~ — m~„ 

<3i x{ 


526 ±34 ±13 


516.0 


mr~ — m~„ 
>R xV 


34.2±3.8=F 0.1 


37.6 



9. Determination SUSY Model 
Parameters 

The ultimate goal is to pin down the SUSY model 
parameters. The SUSY parameter-fitting package Fit- 
tino [13] has been used to perform the calculations 
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of the model parameters from all the endpoints mea- 
sured; m 



edge edge thr 
5 m,,„ , m«o r . 



edge 



edge edge 
'"ig(high)' " 1 tt J 



Wt2(5r)- The fitted mSUGRA parameters can be 
seen m Table [Hi] for the SU3 point. Already with 
1 fb _1 of data m and m 1 / 2 parameters can be de- 
rived reliably, however determination of tan(/3) and 
Aq is more problematic since no information from the 
Higgs sector is available at low luminosity. 



Table III Fitted mSUGRA parameters for the SU3 point 
with 1 fb" 1 at 14 TeV. Sgn(/i) = +. 



Parameter 


SU3 Fitted 


SU3 Truth 


tan(/3) 


7.4 ±4.6 


6 


m (GeV) 


98.5 ±9.3 


100 


m 1/2 (GeV) 


317.7 ±6.9 


300 


A (GeV) 


445 ± 408 


-300 



10. Conclusions 

If SUSY particles exist in nature at sub- TeV mass 
range the ATLAS and CMS experiments will open up 
a discovery window for new physics beyond the Stan- 
dard Model. Search strategies have been developed 
to cover a broad range of signatures expected from 
SUSY models. Understanding the detector perfor- 
mance and controlling the SM backgrounds via data- 
driven methods are important tasks for early SUSY 
analyses. Once a signature consistent with SUSY has 
been established then various measurements will have 
to be combined to reconstruct the SUSY mass spec- 
trum and to understand the SUSY breaking mecha- 
nism. 

Early SUSY searches with 100/200 pb _1 integrated 
luminosity at 10 TeV center of mass energy suggest 
that we are already sensitive to an unexplored area 
in less than a year running. Exciting times are ahead 
with the upcoming start-up of the LHC. 
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